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in : 

O ■ 1. A BRIEF HISTORY 

o ■ 

CN ' The "Falicov-Kimball model" was first considered by Hubbard and Gutzwiller in 1963- 

65 as a simplification of the Hubbard model. Falicov and Kimball introduced in 1969 
a model that included a few extra complications, in order to investigate metal-insulator 
phase transitions in rare-earth materials and transition-metal compounds. Experimental 
. data suggested that this transition is due to the interactions between electrons in two 

electronic states: non-localized states (itinerant electrons), and states that are localized 
■ around the sites corresponding to the metallic ions of the crystal (static electrons). 

A tight-binding approximation leads to a model defined on a lattice (the crystal) and 
two species of particles are considered. The first species consists of spinless quantum 
fermions (we refer to them as "electrons"), and the second species consists of localized 
holes or electrons ("classical particles"). Electrons hop between nearest-neighbor sites 
but classical particles do not. Both species obey Fermi statistics (in particular, the Pauli 
exclusion principle prevents more than one particle of a given species to occupy the same 
site). Interactions are on-site and thus involve particles of different species; they can be 
^ ■ repulsive or attractive. 

The very simplicity of the model allows for a broad range of applications. It was studied 
' i n the context of mixed valence systems, binary alloys, and crystal formation. Adding a 

magnetic field yields the flux phase problem. The Falicov-Kimball model can also be 
viewed as the simplest model where quantum particles interact with classical fields. 

The fifteen years following the introduction of the model saw studies based on approx- 
imate methods, such as Green's function techniques, that gave rise to a lot of confusion. 
A breakthrough occurred in 1986 when Brandt and Schmidt, and Kennedy and Lieb, pro- 
posed the first rigorous results. In particular, Kennedy and Lieb showed in their beautiful 
paper that the electrons create an effective interaction between the classical particles, and 
that a phase transition takes place for any value of the coupling constant, provided the 
temperature is low enough. 

Many studies by mathematical-physicists followed and several results are presented in 
this short survey. Recent years have seen an increase of interest by condensed matter 
physicists. Only a few references are given here, see the "bibliographical notes" below. 
We encourage interested readers to consult the reviews (fH E] . 

2. Mathematical setting 

2.1. Definitions. Let Ac^ denote a finite cubic box. The configuration space for the 
classical particles is 

tt A = {0, 1} A = {u = (uj x ) : x <E A, and u x = 0, l}, 
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where ui x = or 1 denotes the absence or presence of a classical particle at the site x. The 
total number of classical particles is N c (uj) = YlxeA^x- The Hilbert space for the spinless 
quantum particles ("electrons") is the usual fermionic Fock space 

|A| 

^A = (J) T~Ca,n, 

N=0 

where Ha,n is the Hilbert space of square summable, antisymmetric, complex functions 
<I> = . . . , xjv) of N variables x.- L € A. Let a x and a x denote the standard creation and 

annihilation operators of an electron at x; recall that they satisfy the anticommutation 
relations 

{a x ,a y } = 0, {a x ,al} = 0, {a x ,a\} = 5 xy . 

The Hamiltonian for the Falicov-Kimball model is an operator on Ts. that depends on the 
configurations of classical particles. Namely, for uj € Qa, we define 

H a (uj) = - a l a y ~ U Y^ w *4 a z- 

\x-y\=l 

The first term represents the kinetic energy of the electrons. The second term repre- 
sents the on-site attraction (U > 0) or repulsion (U < 0) between electrons and classical 
particles. 

The Falicov-Kimball Hamiltonian can be written with the help of a one-body Hamil- 
tonian h,A, which is an operator on the Hilbert space for a single electron £ 2 (A). Indeed, 
we have 

Ha(uj) = ^ h xy (u)a x a y . 

x,y<=A 

The matrix h\(uj) = (h xy (u)) is the sum of a hopping matrix (adjacency matrix) 
and of a matrix v\(u) that represents an external potential due to the classical particles. 
Namely, we have 

h xy (co^ — t xy Uuj x S xy , 

where t xy is one if x and y are nearest-neighbors, and is zero otherwise. The spectrum of 
t\ lies in (—2d, 2d), and the eigenvalues of v\(u>) are — U (with degeneracy N c (u)) and 
(with degeneracy |A| — N c (u>)). Denoting Xj(A) the eigenvalues of a matrix A, it follows 
from the minimax principle that: 

Xj(A)- \\B\\ ^ \j(A + B) < \j(A) + \\B\\. 

Let \\(oj) ^ A2(cj) ^ ... ^ A|A|(w) be the eigenvalues of h\ (ui). Choosing A = v\(u) 
and B = t\ in the inequality above, we find that for U > 0, 

-U-2d< Xj(u) <-U + 2d for j = 1, ... , N c (lo), 

-2d < Xj(u) < 2d for j = N c (lo) + 1, . . . , |A|. 

In particular, for any configuration uj and any A, 

Spec h A (uj) C (-U - 2d, -U + 2d) U (-2d, 2d). 

Thus for U > Ad the spectrum of h\(u) has the "universal" gap (—U + 2d, —2d). A similar 
property holds for U < —Ad. 
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2.2. Canonical ensemble. A fruitful approach towards understanding the behavior of 
the Falicov-Kimball model is to fix first the configuration of the classical particles and 
to introduce the ground state energy E\(N e ,oj) as the lowest eigenvalue of H\(uj) in the 
subspace TC\ t N c : 

E A (N e ,u;) = inf {*\H A (u) =VA». 
*6Wa,jv c ,||*||=i j^i 

A typical problem is to find the set of ground state configurations, i.e. the set of configu- 
rations that minimize E\(N e ,tu) for given N e and N c = N c (u>). 

In the case U > 4d and A c = N c (u), the ground state energy E\(N c (u>),u>) has a 
convergent expansion in powers of U : 

e a (n^ M = -un^ )+ J:^ zt Y. (- 1 )"'* : ' ,, („ i({ V ) )-i)- (D 

k^2 X!,...,x k eA 7 

\xi— aii_i]=l 
0<m({x.i})<k 

where m(xi, . . . , x^) is the number of sites Xi with u Xi = 0. The last sum also includes 

2d 
U-4d J 



the condition \xk — x\\ = 1. Simple estimates show that the series is less than jj 2 ^ a , N c {uj). 



The lowest order term is a nearest-neighbor interaction, 

{x,y}:\x-y\=l 

that favors pairs with different occupation numbers. Formula Q is the starting point for 
most studies of the phase diagram for large U. A similar expansion holds for U < —Ad 
and A c = |A| - JV c (u>). 

Phase diagrams are better discussed in the limit of infinite volumes where boundary 
effects can be discarded. Let f2 per be the set of configurations on Z d that are periodic in 
all d directions, and O pcr (>o c ) C O per be the set of periodic configurations with density p c . 
For oj £ O per and p c E [0, 1], we introduce the energy per site in the infinite volume limit 
by 

e{p e ,uj)= lim r^E A (N c ,u). (2) 

Here, the limit is taken over any sequence of increasing cubes, and A c = Lpe|A|J is the 
integer part of p e | A. | . Existence of this limit follows from standard arguments. 

In the case of the empty configuration uj x = 0, we get the well-known energy per site of 
free lattice electrons: For k G [— ir, 7r] d , let e(k) = — Ylt=i cos ^v; then 

e(p e ,u> = 0) = — I e{k)dk, 

( 2vr ) Je(k)<e F (p c ) 

where £F(Pe) is the Fermi energy, defined by 

1 



dife. 



(2vr) d Je(k)<e F ( Pc ) 

The other simple situation is the full configuration oj x = 1, whose energy is e(p e ,u> = 1) = 
e(p c ,iv = 0) - Up e . 

Let e(p e , p c ) denote the absolute ground state energy density, namely, 

e(p e ,p c )= inf e(p e ,u). 

Notice that e(p e ,u>) is convex in p c , and that e(p e ,p c ) is the convex envelope of {e(p e ,uj) : 
co G f2 per (p c )}. It may be locally linear around some (p e ,p c ). This is the case if the 
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infimum is not realized by a periodic configuration. The non-periodic ground states can 
be expressed as linear combinations of two or more periodic ground states ("mixtures"). 
That is, for 1 ^ i ^ n there are cq ^ with Oj = 1, G $7 per , and pi' , such that 

i i 

and 

i 

The simplest mixture is the "segregated state" for densities p c < p c : Take u/ 1 ) to be 
the empty configuration, uj^ to be the full configuration, p^ = 0, p^ = and 02 = 

1 - «1 = Pe- 
lf d ^ 2, a mixture between configurations cjW can be realized as follows. First, 

partition 7L d into domains D\ U • • • U D n such that ^4 — > Qj and — > as A f r L d . 

(i) 

Then define a non-periodic configuration u by setting uj x = ujx for x £ D{. See the 
illustration in Fig. ^ The canonical energy can be computed from (J2J), and it is equal to 

n 

e(p e ,w)= inf Vaie(pW,wW). 



=Pe j = l 



Furthermore, the infimum is realized by densities p\ such that there exists p e with 
p e (p e ,io^) = pe for all i (See Q below for the definition of p c {pe, 
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Figure 1. A two-dimensional mixed configuration formed by periodic con- 
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figurations of densities 0, i, and 1 



We define the canonical ground state phase diagram as the set of ground states u> (either 
a periodic configuration or a mixture) that minimize the ground state energy for given 
densities p e , p c : 

G ca , n (p e ,Pc) = : e(p e ,uj) = e(p c ,p c ) and p c (u) = p c }. 

2.3. Grand-canonical ensemble. Properties of the system at finite temperatures are 
usually investigated within the grand-canonical formalism. The equilibrium state is char- 
acterized by an inverse temperature (3 = l/k^T, and by chemical potentials p e , p c , for the 
electrons and for the classical particles respectively. In this formalism the thermodynamic 
properties are derived from the partition functions 

Z A (P,p e ,u) = Tr^ A e-^ A H-^iVA] , 

Z A ((3,p e ,p c )= e^ N ^ Z A ((3,p e ,u). ( 3 ) 
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Here, N A = YlxeA ^a x is the operator for the total number of electrons. We then define 
the free energy by 

F A (f3,p e ,p c ) = - ± log Z A ((3,p e ,p c ). 

The first partition function in (J3J) allows to introduce an effective interaction for the 
classical particles, mediated by the electrons, by 

F A ((3,p e ,p c ,uj) = -p c N c (u) - j$ log Z A (P,p e ,uS). 

It depends on the inverse temperature (3. Taking the limit of zero temperature gives the 
corresponding ground state energy of the electrons in the classical configuration lo, 

E A (p e ,p c ,Cj) = )hlL F A (f3,p e ,pl c ,u) = -p c N c (u) + V] (Xj(u)-p e ). 
p^oo ' — * 

j:\ J (uj)<fi c 

Notice that F A and E A are strictly decreasing and concave in p e ,p c (E A is actually linear 
in p c ). We also define the energy density in the infinite volume limit by considering a 
sequence of increasing cubes. For to G f2 per , 

e(p e ,p c ,u;) = lim ^E A (p, e ,p c ,u)). 

The corresponding electronic density is 

p e (p e ,u) =^lim d #{j : Xj(uj) < p c } = -^j-e{p e ,p c ,uj), (4) 

and the density of classical particles is p c {u) = lim A N c (lu)/\A\. One can check that 
canonical and grand-canonical energies are related by 

e(/x e , ix c ,uj) = e(p e ((j, e ,u),u) - p c p e (p Cl io) - p c p c (u). (5) 

Given (p e ,p c ), the ground state energy density e(p e ,p c ) is defined by 

e(p e ,p c )= inf e(p e ,p, c ,u). 

The set of periodic ground state configurations for given chemical potentials p e , p c is the 
grand- canonical ground state phase diagram: 

G gc (pe,p c ) = {uj E f2 pcr : e(p e ,p c ,uj) = e(p c ,p c )}. 

It may happen that no periodic configuration minimizes e(p c , p c , uj) and that G gc (p e , p c ) = 
0. Results suggest that G gc (p e , p c ) is nonempty for almost all p e , Pc, however. 

The situation simplifies for U > 4d and p e € (-U + 2d, —2d). Since p e belongs to the 
gap of h A (uj), we have p e (p e ,u>) = p c (w), and 

e(p c ,p c ,uj) = e(p c (uj),uj) - (p c + p c )p c {u). 

Thus G gc (p e , p c ) is invariant along the line p e + p c = const (for p e in the gap). 

2.4. Symmetries of the model. The Hamiltonian H A clearly has the symmetries of the 
lattice (for a box with periodic boundary conditions, there is invariance under translations, 
rotations by 90°, and reflections through an axis). More important, it also possesses 
particle-hole symmetries and these are useful since they allow to restrict investigations to 
positive U and to certain domains of densities or chemical potentials (see below). 

• The classical particle-hole transformation tu x i— ► ZJ X = 1 — lv x results in 

HV(oJ) = H a u (oj) - UN A , 



6 CHRISTIAN GRUBER, DANIEL UELTSCHI 

and N c (uJ) = |A| - N c (u). It follows that E^(N e ,u) = E~ u (N c ,oj) - UN C , and 

Gcan(Pc,Pc) = {cJ : OJ G G^ an (p e , 1 - /5 C )}, 
G~^(p e ,p c ) = {iO :UJ e Gg C (// e ~ ?/, ~fJ>c)}- 

• An electron-hole transformation can be defined via the unitary transformation 
a x l— * £x&x and a x i— > e^a^, where is equal to 1 on a sublattice, and to —1 on 
the other sublattice. Then 

H%(u)~HZ u (u)-UN c (u), 

and N A i-» |A|-AT a . It follows that ^(|A| - JV e , w) = Ej? (N e ,u) -UN c (w), and 

Gcan(Pe,Pc) = Gca n (l — p e ,Pc), 
Ggc {pe, pc) = Gg C ( — p e , p c — U). 

• Finally, the particle-hole transformation for both the classical particles and the 
electrons give 

H% (57) i-» H a (lo) + IW A + J7iV c (a;) - f/|A|. 
It follows that E%(\A\ - N e ,UJ) = E^(N e ,u) + U(N e + N c (lo) - |A|), and 

G^(pe,Pc) = {uJ-.UJ £ Gcijl ~ pe,l~ Pc)} , 
G^ c (p e ,p c ) = {uJ : id G Gg C (-p e - U,-p c - U)}. 

Any of the first two symmetries allow to choose the sign of U. We assume from now 
on that U ^ 0. The third symmetry indicates that the phase diagrams have a point of 
central symmetry, given by p e = p c = i in the canonical ensemble and p e = p c = — ^ 
in the grand-canonical ensemble. Consequently, it is enough to study densities satisfying 
p c ^ | and chemical potentials satisfying p e ^ — y. 

These symmetries have also useful consequences at positive temperatures. In particular, 
both species of particles have average density | at p e = p c = — ^, for all /?. 

3. The ground state — Arbitrary dimensions 

3.1. The segregated state. What follows is best understood in the limit U — > oo and 
when p c < p c . In this case the electrons become localized in the domain T> A (u>) = {x G A : 
w x = 1} and their energy per site is that of the full configuration, e(p, u = 1) (see Section 
12.2(1 . where p = p e / p c is the effective electronic density. The presence of a boundary for 
V a (lo) raises the energy and the correction is roughly proportional to 

B A (u) = #{(x,y) : x G V A (co) and y G Z d \V A (u)}. 

Theorem 1. 

(a) Let A C 7L d be a finite box, and U > Ad. Then for all lo G fl A , an d a ^ N e sj 
N c (u) = N c , we have the following upper and lower bounds: 

^\e(^)\B A (u) > E A (N e ,u)-N c e(f c ^ = l) > [a(f ) - -y(U)] B A {u). 

Here, a{p) = a(l — p) is strictly positive for < p < 1. 7(17) behaves as ^ for 
large U , in the sense that U^{U) — > 8d 2 as U — > oo. 

(b) For any p e ^ p c that differ from zero, the segregated state is the unique ground 

■ Pc ' 



state if a(Pf-) > j(U), i.e. if U is large enough. 
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The proof of (a) is rather lengthy and we only show here that it implies (b). Let 
b(uj) = liniA ? an£ i notice that b(uj) = for the empty, the full, and the segregated 
configurations; < b(uj) < d for all other periodic configurations or mixtures. Recall that 
p c e(^,uj = 1) is the energy density of the segregated state. For all densities such that 
a(^) > "y(U), and all configurations such that p c {ui) = p c , we have 

e(p e ,u) > p c e(^,uj = 1), 

and the inequality is strict for any periodic configuration. This shows that the segregated 
configuration is the unique ground state. 

3.2. General properties of the grand-canonical phase diagram. We have already 
seen that the grand-canonical phase diagram is symmetric with respect to (— y). 
Other properties follow from concavity of e(p e ,p c )- 

Let uj G G gc (p e ,p c ) \ G gc (p' c ,p' c ) and uj' G G gc (p! e ,p' c ) \ G gc (p e ,p c )- Then 

(i) p e = p' e and p' c > p c imply p c {uj') > p c (w)) 

(ii) p c = p! c and p' c > p e imply p e (pf e ,uj') > p e (p, e ,u), and uj cannot be obtained by 
adding some classical particles to the configuration <J . 

It follows from (ii) that if u = 1 G G gc (p e , p c ), then uj = 1 G G gc (p' e , p! c ) for all p e //g, 
*5 /^c- A similar property holds for the empty configuration. To establish these 
properties, we can start from 

e(p e ,p c ,uj') - e(p e ,p c ,uj) > > e{p' e ,p! c ,J) - e(p c ,p' c ,u;). (6) 

Since e(p e , p c , lo) is concave with respect to p e and linear with respect to p c , we have 

e(p' e ,p' c ,uj) ^ e(p c ,p c ,u;) + {p c - p' c )p e (p e ,uj) + (p c - p' c )p c {uj). (7) 

Using this inequality for both terms of the right side of Q, we obtain the inequality 

(p' e - p e ) [pe{p' e , Uj') - p e {p e ,Uj)] + (p' c - p c ) [p c (uj') - p c {uj)] ^ 0, 

which proves (i) and the first part of (ii). The second part of (ii) follows from 

/He 
p c (p,uj)dp - p c pc(uj). 
-oo 

Indeed, the minimax principle implies that eigenvalues Xj(u>) are decreasing with respect 
to lu (if U ^ 0), so that p c (p c ,uj) is increasing (with respect to uj). Then for any oj" > uj 
and p' e > p e , 

e(p' e ,pc,uj") - e(p'c,Pc,uj) > e(p e ,Pc,uj") - e{p e ,Pc,uj), 

and uj" £ G gc (p e ,p c ) implies uj" G gc (p' c ,Pc)- 

Next we discuss domains in the plane of chemical potentials where the empty, full, and 
chessboard configurations have minimum energy. One easily sees that uj = 1 is the unique 
ground state configuration if p c > 0, or if p c > 2d and p c > —U. Similarly, uj = is the 
unique ground state if p c < —U, or if p e < —U — 2d and p c < 0. For U > 4d, it follows from 
the expansion (0) that the full configuration is also ground state if — U + 2d < p e < —2d 
and pe + Pc + U > , . These domains can be rigorously extended using energy estimates 
that involve correlation functions of classical particles. The results are illustrated in Figs 
|U(C7 < id) and[2([/ > Ad). 

Finally, canonical and grand-canonical phase diagrams are related by the following 
properties: 

(iii) If uj G G gc (p e ,Pc), then uj G G ca , n (p c (p c ,uj), p c (uj)). 
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Figure 2. Grand-canonical ground state phase diagram for U < 4d. Do- 
mains for the empty, chessboard, and full configurations, are denoted in 
light gray, black, and dark gray respectively. 
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Figure 3. Grand-canonical ground state phase diagram for U > 4d. Do- 
mains for the empty, chessboard, and full configurations, are denoted in 
light gray, black, and dark gray respectively. 

(iv) More generally, suppose that u^',...,u^ € G gc (p e , Pc), and consider a mix- 
ture with coefficients a±, . . . ,a n . The mixture belongs to G can (p e > Pc), with p c = 
Y, i acip e (fJ>e,u^) and p c = J2i a iPc(^)- 

To establish (iii), observe that any u' satisfies e(p e , p c , u') ^ e(p e , p c , uj) if cj G G gc (p e , p, c ). 

Let p e = p e (p e ,uj) and p c = p c (u), and let p! e such that p e (p' e ,oj') = p e . By Eqs (jSJ) and 

0, 

e(p e (fi' e ,u'),u') - p c p c (p' e ,uj') - fi c p c (u') ^ e(p e (p e ,u),u) - p e p e (p e ,oj) - p c p c {u). 

Then e(p e ,u') ^ e(p e ,u) for any configuration u/ such that p c {u') = p c . Property (iv) 
follows from (iii) by a limiting argument, because a mixture can be approximated by a 
sequence of periodic configurations. 

Next we describe further properties of the phase diagrams that are specific to dimensions 
1 and 2. 

4. Ground state configurations — Dimension one 

A large number of investigations, either analytical or numerical, have been devoted to 
the study of the ground state configurations in one dimension. One-dimensional results 
also serve as guide to higher dimensions. Recall that symmetries allow to restrict to 
U ^ and p c ^ 5. 
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Most ground state configurations that appear in the canonical phase diagram seem to 
be given by an intriguing formula, which we now describe. Let p e = | with p relatively 
prime to q. Then corresponding periodic ground state configurations have period q and 
density p c = 7 - (r is an integer). The occupied sites in the cell {0, 1, . . . , q — 1} are given 
by the solutions fco, . . . , fc r -l of 

(pkj) = j mod q, ^ j ^ r — 1. (8) 

Note that the first classical particle is located at k$ = 0, and k$, . . . , k p _\ are not in 
increasing order. In order to discuss the solutions of ©, we introduce I = |_pj (the 
integer part of and we write 

q=(t+l)p-s, (9) 

where 1 ^ s ^ p — 1, and s is relatively prime to p. Next, let L(x) denote the distance 
between the particle at x and the one immediately preceding it (to the left). 
Let us observe that if p c = p c , i.e. if r = p, then 

(a) L(kj) = £ for ^ j ^ s — 1 and k j — t = kj +p ^ s . 

(b) L(kj) = i + 1 for s ^ j ^ p - 1 and kj - (£ + 1) = fy_ s . 

Indeed, for pkj = j + nq, Eq. Q implies 

p(kj — I) = j + (n — l)q + (p — s) = j + p — s mod q, 

and 

p(kj — t — 1) = j — s mod q. 

Therefore kj — £ is solution of @ if j + p — s ^ p — 1, while fcj — (■£ + 1) is solution of © 
ifj-s ^ 0. 

These properties show that the configuration defined by (jSJ) is such that L(x) £ {^, ^+1} 
for all occupied x. A periodic configuration such that all distances between consecutive 
particles are either i or £ + 1 is called homogeneous. Let w be a homogeneous configuration 
with period q and density p c = ^, and let xo < •■■ < x p ^\ be the occupied sites in 
{0, 1, . . . , q — 1}. We introduce the derivative cu' of lo as the periodic configuration with 
period r defined by (see Fig. 0J) 

, fl if = 4 
Wi [0 if L(aJi) = f + 1. 

A configuration is most homogeneous if it can be "differentiated" repeatedly until the 
empty or the full configuration is obtained. 
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Figure 4. The configuration oj given by the formula Q with q = 24 and 
p = 7, and its derivative a/. Notice that £ = 3 and s = 4. 

Let w be the homogeneous configuration from (jHJ) and a/ be its derivative. Using the 
same arguments as for properties (a) and (b) above, and the fact that s is relatively prime 
to p, we obtain: 
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(c) Let k' , . . . , k' p _ 1 be the solutions of 

(skj) = j mod p. 

Then (k' , . . . , k' p _i) is a permutation of (0, 1, ... ,p — 1). Further, kj — 1 
for ^ j ^ s — 1, and kj — 1 = for s ^ j ^ p — 1. 

Consider the periodic configuration with period p where sites k' Q , . . . , k' s _ 1 are occupied 
and sites k' s , . . . , k' p _ 1 are empty. Since k' = 0, this configuration is precisely the derivative 
u>' of lu. Iterating, these properties prove that the solutions of (JSJ) are most homogeneous. 

One of the most important result in one dimension is that only most homogeneous 
configurations are present in the canonical phase diagram, for U large enough and for 
equal densities p e = p c . 



-. There exists a constant c such that for U > cA q , 

g •* 



Theorem 2. Suppose that p e = p t 

the only ground state configuration is the most homogeneous configuration, given by 
(together with translations and reflections). 

This theorem was established using the expansion (^Q) of E^(N e , to) in powers of U . It 
suggests a devil's staircase structure with infinitely many domains. However, the number 
of domains for fixed U could still be finite. Results from Theorem [2] are illustrated in Fig. 

El 







-t/ / 






LO = 1 


LO = ^ 




-u- 









Figure 5. Grand-canonical ground state phase diagram in one dimension 
for U > 4 and p, e in the universal gap. Chessboard configurations occur in 
the black domain. Dark gray oblique domains correspond to densities g, 
4' 1> 5' I' I> f - Total width of these domains is of order f/ -1 . 



For small U, on the other hand, one can use a (non-rigorous) Wigner-Brillouin degen- 
erate perturbation theory (a standard tool in band theory). Let p e = | with p relatively 
prime to q, and to be a periodic configuration with period nq, n € N. Then for U small 
enough (U ^ -), we obtain the following expansion for the ground state energy 

e(p e ,u) = -f sinvrpe - Up e p c {u) - ^^-U 2 \ log U\ + 0(U 2 ), (10) 
where lv(p e ) is the "structure factor" of the periodic configuration u, namely 
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This expansion suggests that the ground state configuration can be found by maximizing 
the structure factor. The following theorem holds independently of U. 

Theorem 3. Let p e = |. There exist r% ^ | and r2 ^ ^ such that the configurations 
maximizing the structure factor are given as follows: 

(a) For p c = ^ with r\ ^ r ^ r2, use the formula (JHJ). 

(b) For p c £ (£ u>ii/i ri ^ r ^ r2 — 1, i/ie configuration is a mixture of those 

f or Pc= r - q and p c = £±i . 

(c) For /3 C E (0, — ) ; £/ie configurations are mixtures of lo = and i/iai /or p c = —. 
For p c E (£2-, 1) ; £/ie configurations are mixtures of uj = 1 and i/iai /or /5 C = y . 

Some insight for low densities is provided by computing the energy of just one classical 
particle and one electron on the infinite line, and to compare it with two consecutive 
classical particles and two electrons. It turns out that the former is more favorable than 

the latter for U > -7= ~ 1.15, while "molecules" of two particles are forming when U < 

V3 

-Jg. Smaller U shows even bigger molecules for p c = np e , and n-molecules are most 
homogeneously distributed according to the formula (JHJ). It should be stressed that the 
canonical ground state cannot be periodic if U is small and p c ^ [|, |], which is different 
from the case of large U. 

Only numerical results are available for intermediate U. They suggest that configura- 
tions occuring in the phase diagram are essentially given by Theorem |21 (together with 
the segregated configuration) . This is sketched in Fig. H2 where bold coexistence lines for 
p c > —U — 2 and p e < 2 represent segregated states. 




Figure 6. Grand-canonical ground state phase diagram for U ~ 0.4. En- 
larged are domains for p e = I and | , with the same densities p c = I , I , % ■ 



5. Ground state configurations — Dimension two 

We discuss the canonical ensemble only, but many results extend to the grand-canonical 
ensemble. Recall that G can (^, \) consists of the two chessboard configurations for any 
U > 0, and that segregation takes place when p e ^ p c , providing U is large enough 
(Theorem Other results deal with the case of equal densities, and for U large enough. 

Theorem 4. Let p c = p c = p ^ \- 

(a) If p G {7;, I , |, j, |, 4, ^j-, |}, then for U large enough, the ground state configura- 
tions are those displayed in Fig. Q // p = with integer n, then for U 
large enough (depending on p), the ground state configurations are periodic. 
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(b) If p is a rational number between g and |, then for U large enough (depending on 
the denominator of p), the ground state configurations are periodic. Further, the 
restriction to any horizontal line is a one- dimensional periodic configuration given 
by © , and the configuration is constant in either the direction ( } ) or ( _} 1 ) . 

(c) Suppose that U is large enough. If p € (g>fy)> the ground state configurations 
are mixtures of the configurations p = \ and p = of Fig. If p G (g, |), the 
ground state configurations are mixtures of the configurations p = \ o,nd p = | . If 
p G (|, j), the ground state configurations are mixtures of the configurations p = | 
and p = \ - 

p c is presented in Fig. |HJ 



The canonical phase diagram for p c 
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Figure 7. Ground states configurations for several densities. Occupied 
sites are denoted by black circles, empty sites by white circles. Lines are 
present only to clarify the patterns. 

The situation for densities p ^ \ that are not mentioned in Theorem |1] is unknown. 
All these periodic configurations are present in the grand-canonical phase diagram as well. 
Item (b) suggests that the two-dimensional situation is similar to the one-dimensional one 
where a devil's staircase structure may occur. Let us stress that no periodic configurations 
occur for large U and densities p e = p c in the intervals (g, ^-), (^,|), and (§,7). This 
resembles the one-dimensional situation, but for small U. 
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Figure 8. Canonical ground state phase diagram in two dimensions for 
U >8. 



Bibliographical notes 

The Falicov-Kimball model was introduced in [5] , and the first rigorous results appeared 
in and ^U]. A simple derivation of the expansion (^Q) using Cauchy formula can be found 
in [7]. It can be extended to positive temperatures with the help of Lie-Schwinger series 
[2]. Segregation in arbitrary dimension (Theorem ^) was proposed in 5 . Domains in 
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the grand-canonical phase diagram for the empty, full, and chessboard configurations, 
were obtained using lower bounds involving correlation functions of classical particles, see 
[7] for references. Eq. ()1U|) for the ground state energy for U small was derived using 
Wigner-Brillouin degenerate perturbation theory in [I] . For the results in two dimensions 
(Theorem we refer to [5] and references therein. Results on the ground state for 
jj, e in the universal gap have been extended to positive temperatures using "quantum 
Pirogov-Sinai theory" (see the review in this Encyclopedia). Further studies and extensions 
of the Falicov-Kimball model include interfaces, non-bipartite lattices, bosonic particles, 
continuous fields instead of classical particles, magnetic fields, correlated and extended 
hoppings. An extensive list of references can be found in the reviews jSl|7||5]. 

See also: Equilibrium Statistical Mechanics. Quantum Statistical Mechanics. Fermionic 
Systems. Hubbard Model. Pirogov-Sinai Theory. 
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